Orbital-controlled magnetic transition between gapful and gapless phases 
in the Haldane system with t2g-orbital degeneracy 
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In order to clarify a key role of orbital degree of freedom in the spin 5=1 Haldane system, we 
investigate ground-state properties of the t2g-orbital degenerate Hubbard model on the linear chain 
by using numerical techniques. Increasing the Hund's rule coupling in multi-orbital systems, in 
general, there occurs a transition from an antiferromagnetic to a ferromagnetic phase. We find 
that the antiferromagnetic phase is described as the Haldane system with spin gap, while in the 
ferromagnetic phase, there exists the gapless excitation with respect to orbital degree of freedom. 
Possible relevance of the present results to actual systems is also discussed. 
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Since the Haldane conjecture,^ quasi-one-dimensional 
Heisenberg antiferromagnets with spin 5>l/2 have at- 
tracted considerable attention in the research field of con- 
densed matter physics. At first, it was unacceptable that 
the integer-spin chain exhibits an energy gap in the spin 
excitation and the spin correlation decays exponentially 
in contrast to the half-odd-integer-spin chain. However, 
the prediction has been confirmed to be correct after 
examined by theoretical, ^'"^'^ numericalj^*^ and experi- 
mental investigationsiSiSiiii Especially, the valence-bond- 
solid (VBS) picture has clarified the microscopic mecha- 
nism of the gapped ground state, ^ and the Haldane gap 
has been actually observed for spin 5*=! compounds such 
as Ni(C2HgN)2N02(C104) abbreviated NENP.^ It has 
been also found that when a magnetic field is applied 
to NENP, finite magnetic moment begins to grow at a 
critical field corresponding to the Haldane gapiS 

Now the existence of the Haldane gap is confirmed 
in the spin S=l system, but in LiVGe2 06, a quasi-one- 
dimensional spin 5'=1 system, antiferromagnetic (AFM) 
transition occurs at a Neel temperature Tn=22K and the 
expected Haldane gap is absent or strongly suppressed. 
To understand the suppression contrary to the Hal- 
dane conjecture, a scenario has been proposed based on 
biquadratic interaction of neighboring spins. In fact, 
the susceptibility has been reproduced by the bilinear- 
biquadratic S=\ Heisenberg chain, in combination with 
the effect of impurity contribution^ The biquadratic in- 
teraction has been derived by the fourth-order perturba- 
tion in terms of the electron hopping,-^ assuming a sig- 
nificant value of the level splitting in orbitals of V'^+ 
ion. However, the splitting has been found to be much 
smaller than previously considered,^'* indicating that or- 
bital fiuctuations may play a crucial role in LiVGe2 06. 

Here let us reconsider the Haldane system from an elec- 
tronic viewpoint. As recognized in the VBS picture, spin 
5=1 is decomposed into a couple of 5=1/2 electrons. 
In actual systems, the spin 5=1 state is stabilized by 
the Hund's rule coupling among electrons in different or- 
bitals. ForNi2+ ion in NENP, two electrons in eg orbitals 
form 5=1 and orbital degree of freedom is inactive. On 
the other hand, ion contains two electrons in three 



i2g orbitals. When the level sphtting among orbitals 
is small, as actually observed in LiVGe206, orbital degree 
of freedom remains active after forming 5=1, in contrast 
to Ni^+ ion. Significance of t2g-orbital degree of free- 
dom has been also emphasized to understand the peculiar 
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situation may occur in one-dimensional ruthenium com- 
pounds, since the low-spin state of Ru'*+ ion with Adt^ 
configuration includes two holes. In the Haldane system, 
we have learned that distinctive features originate from 
strong quantum effects of the spin. We believe that it is 
fascinating to consider further orbital degree of freedom 
in the Haldane system. In particular, it is important to 
clarify how orbital degeneracy affects magnetic proper- 
ties such as the gapful nature of the spin-only system. 

In this paper, we investigate ground-state properties 
of the <2g-orbital degenerate Hubbard model on the lin- 
ear chain by numerical techniques such as the Lanczos 
diagonalization and the density matrix renormalization 
group (DMRG) methodiiS, When the Hund's rule cou- 
pling is small, an AFM phase with ferro-orbital (FO) or- 
dering appears. This phase is well described by the spin 
5=1 AFM Heisenberg model and thus, the Haldane gap 
exists. On the other hand, increasing the Hund's rule 
coupling, there occurs a ferromagnetic (FM) phase with 
antiferro-orbital (AFO) correlation, in which the system 
is described by the pseudo-spin T=l/2 AFO Heisenberg 
model with the gapless orbital excitation. Namely, the 
low-energy physics drastically changes between gapful 
and gapless phases as a result of orbital ordering. 

Let us consider a chain of ions including four electrons 
per site, where the local spin 5=1 state is expected to 
form in i2g orbitals. Here we envisage Ru^^ chain, but 
due to the electron-hole symmetry, the results are also 
applicable to chain with two electrons in orbitals 
per site. Note that the direction of the chain is set to 
be the x-axis, but the result does not depend on the 
chain direction due to the cubic symmetry. The orbital 
degenerate Hubbard model is given by 



H = 
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where di-yo- is an annihilation operator for an electron 
with spin a in the orbital 7 {=xy, yz, zx) at site i, t^-^i is 
the nearest-neighbor hopping between adjacent 7 and 7' 
orbitals, and Pi-ya=dl^^dija- Note that hopping ampli- 
tudes are given as txy^xy=tzx,zx=t and zero for other com- 
binations of orbitals. Since we consider the chain along 
the a;-axis, there is no hopping for yz orbitals. Hereafter, 
t is taken as the energy unit. In the interaction terms, 
U is the intra-orbital Coulomb interaction, U' the inter- 
orbital Coulomb interaction, J the inter-orbital exchange 
interaction (the Hund's rule coupling), and J' is the pair- 
hopping amplitude between different orbitals. Note the 
relation U—U'+J+J' due to the rotation symmetry in 
the orbital space and J'=J is also assumed due to the 
reality of the wavefunction. 

First we show the results of the Lanczos diagonaliza- 
tion for the 4-site chain with the anti-periodic boundary 
condition. Although the system size is very small and it 
is necessary to pay due attention to a relation with actual 
materials, there are advantages that we can quickly accu- 
mulate lots of results and easily gain insight into ground- 
state properties. In Fig. the ground-state phase dia- 
gram is shown in the {Ucs, J) plane with Ucs—U'^J. 
The phase boundary is determined by comparing the en- 
ergies for the phases of S^^^=0^4:, where S^^^ is the z- 
component of total spin. Here two types of the ground 
state are observed. In the region of small J, the lowest- 
energy state is characterized by S^^^=0, indicating the 
spin singlet state regarded as the one-dimensional AFM 
phase. On the other hand, the maximally spin-polarized 
FM phase appears in the region of large J, where the 
energies for the phases with 5*^0^=0^4 are degenerate. 

In order to obtain an intuitive explanation for the ap- 
pearance of both phases, we show the electron config- 
uration in t2g orbitals for the AFM and FM phases in 
the insets of Fig. 1. In the AFM phase, since there is 
no hopping between adjacent yz orbitals along the a:- 
axis, xy and zx orbitals should be singly occupied to 
gain the kinetic energy by electron hopping, while yz or- 
bital is doubly occupied and localized. The Hund's rule 
coupling stabilizes the local spin S—1 state consisting of 
two electrons with parallel spin in xy and zx orbitals, 
which are antiferromagnetically coupled between nearest 
neighbor sites. Note that the electron configuration in 
the AFM phase is well understood from the VBS pic- 
ture: Due to the hopping connection, it is natural that 
a spin singlet pair should be formed on adjacent xy and 
zx orbitals. Considering that the spin singlet pair is al- 
ternately assigned to xy or zx orbital, the uniform VBS 
state is constructed. On the other hand, in the FM phase, 
three up-spin electrons occupy different orbitals due to 
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FIG. 1: Ground-state phase diagram for the 4-site chain. 
Insets indicate schematic views for electron configurations in 
the AFM and FM phases. Oval enclosures represent the spin 
singlet pair for the VBS picture. 

the Pauli principle. One down-spin electron can be ac- 
commodated in any of three orbitals in principle, but in 
order to gain the kinetic energy, down-spin electrons oc- 
cupy xy or zx orbital alternately. By introducing 
pseudo-spin operators representing xy and zx orbitals, 
the AFO structure is expressed as an alternating corre- 
lation of the pseudo-spin. 

From the above discussion on the electron configura- 
tion, it is useful to consider the effective model in the 
strong-coupling region to clarify ground-state properties. 
For the AFM phase, the effective model is given by 

-f^AFM = J AFM ^(Sj • Sj+i — 1), (2) 

i 

where Si indicates the S—1 operator and the AFM cou- 
pling Jafm is expressed as JAFM=2i^/(J7 -I- J). On the 
other hand, the effective model in the FM phase is given 
by the T=l/2 AFO Heisenberg model 

Hfm = 2Jfm ^(T, • T,+i - 1/4), (3) 

i 

with JpM—2t'^/Ucff- We define pseudo-spin operator 
as Tj=(l/2) Y.j,-y',cr dlja^u'd'iY'^^ where cr=(cr^, cTj^, a^) 
are Pauli matrices and the summation for 7 and 7' is 
taken for xy and zx orbitals. Here we arrive at the fol- 
lowing interesting conclusion: The AFM phase of the 
t2g-orbital degenerate Hubbard model is the gapped Hal- 
dane phase, well described by the spin S'=l AFM Heisen- 
berg model i?AFM- On the other hand, the FM phase 
is described by the pseudo-spin T=l/2 AFO Heisenberg 
model HpM with the gapless excitation. 

We believe that the characteristics grasped within the 
small 4-site chain are physically meaningful, but it is 
highly required to confirm them without finite size ef- 
fects. In order to consider numerically the problem in 
the thermodynamic limit, we employ the infinite-system 




FIG. 2: Ground-state energy per site as a function of Uctt 
for J=6 in the thermodynamic limit. Note that the origin 
of the energy is shifted (see the maintext). Dashed and solid 
curves denote exact ground-state energies for the S=l AFM 
and T=l/2 AFO Heisenberg chains, respectively. 



DMRG method with open boundary conditionii& Note 
that one site consists of three t2g orbitals and the num- 
ber of bases is 64 for one site, indicating that the matrix 
size becomes very large. In the present calculation, the 
number of states kept for each block m is up to r7i=160 
and the truncation error is estimated as 10""" at most. 

In Fig. 121 the ground-state energy is plotted for S^^^^=0 
and N (the number of sites, i.e., the maximally polarized 
state). Here we fix J— 6 to consider the strong-coupling 
region. Note that the energy is shifted by the constant 
of on-site Coulomb interactions 5U' + U — 3J for com- 
parison with the effective model. In the AFM phase, 
the lowest-energy state is characterized by S^^^=0 (solid 
circle), which is in agreement with the exact energy of 
Hafm (dashed curve). On the other hand, the ener- 
gies for S^^^=0 and N are degenerated in the FM phase. 
The DMRG results (open circle) are compared with the 
exact energy of i?FM (solid curve) The agreement is 
fairly well, except for a little discrepancy for small UcS, 
since the strong-coupling limit is appropriate for large 
Ucs- Therefore, we conclude that the ground-state en- 
ergy is well reproduced by the effective model in the 
strong-coupling region. Note that the phase transition 
takes place at a critical point C/c— 33, which is similar to 
the value of the 4-site chain C/c(-^=4)~35. 

In order to verify that ground-state properties are de- 
scribed by the effective model from the microscopic view- 
point, it is quite important to investigate correlation 
functions. For the AFM state, we evaluate the spin corre- 
lation function, CgpinC?, fc) = (5|5f }, the string correlation 

function,!^ Cstring(j, k) = {SI exp[i7r X^fj/ and 
the orbital correlation function, Coihitai{j, k)={TjT^) . 
Here (■ • ■) indicates the average using the ground-state 
wavefunction. In Fig. |3Ja), we show correlation func- 
tions for ;7off =60. We also depict Cspin of Hafm obtained 
by the quantum Monte Carlo (QMC) method with the 
loop algorithm^ We observe that Cspin rapidly decays 
in accordance with the QMC results and the string long- 
range order exists, indicating that the characteristics of 
Hafm sue reproduced and the system has the gapful na- 
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FIG. 3: (a) Spin, string, and orbital correlation functions 
measured from the center of the chain for the AFM state. 
Here r denotes the distance from the center position. Solid 
diamonds denote the spin correlation function of -ffAPM. (b) 
Correlation functions for [/eff=34. (c) The linear-log plot of 
the spin correlation function in Fig. Ola) and the two-orbital 
model, (d) The log-log plot of the orbital correlation function 
for the FM state. 



ture with respect to the spin. Note that the orbital has no 
effect on the low-lying excitation due to the localization 
of yz orbital, as shown by dashed line in Fig. Ufa). For 
C^cff=34 (near the phase boundary), the behavior of cor- 
relation functions in the bulk is also described by Hafm , 
as shown in Fig. Elb). Note that anomaly due to open 
boundary condition appears just at the chain edge, since 
an orbital singlet dimer state with FM spin structure is 
stabilized to gain energy locally at the chain edge. 

In Fig.|3fc), we show Cgpin in the linear- log scale, indi- 
cating evidence of the exponential decay, although there 
is slight deviation from that of the QMC result for Hafm ■ 
In the DMRG calculation for H, the truncation error is 
about 10^^ for to=160 and the energy is evaluated in this 
accuracy. However, in order to evaluate Cspin in the same 
accuracy, it is necessary to increase m further and check 
the convergence with respect to m. Unfortunately, it is 
quite difficult to perform such tasks due to the large ma- 
trix size. Instead, here we introduce a two-orbital Hub- 
bard model composed of xy and zx orbitals, in order to 
reduce the matrix dimension effectively. When we com- 
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pare the results of the effective two-orbital model and the 
original Hamiltonian for m—80 and 160, good agreement 
is obtained as shown in Fig.Ofc), indicating that the two- 
orbital model well reproduces ground-state properties of 
H in the AFM region. Further increasing m up to m=300 
in the two-orbital model to improve the precision of the 
finite-system DMRG method, we confirm that Cspin of 
the two-orbital model converges to that of -ff afm • Then, 
it is highly believed that the same convergence should be 
observed for the original t2g-orbital model. 

For the FM state, we evaluate Corbitai- Note that Cgpin 
is constant depending on S^^^ in the maximally spin- 
polarized FM phase. In Fig. I^d), we observe the power- 
law decay of Corbitai, which is characteristic of the half- 
odd-integer-spin chain. Namely, the system is described 
by ^^FM and quantum fluctuation with respect to the or- 
bital gives the gapless excitation. Further, the spin-wave 
excitation presents due to FM ordering and thus, both 
spin- and orbital- wave excitations are expected to exist. 

Now we discuss relevance of our results to actual sys- 
tems. We have performed the DMRG calculations for 
large value of J=6, to discuss the behavior in the strong- 
coupling region. However, we have confirmed that for 
the small 4-site chain, spin and orbital structures show 
characteristic behavior in each phase even in the weak- 
coupling region with small UcS and J. In addition, as 
mentioned above, the phase boundary for J=6 is almost 
the same as that obtained in the 4-site chain. Namely, 
it is expected that the phase diagram in Fig. ^ does not 
change so much even in the thermodynamic limit, al- 
though the AFM region may somewhat extend. If we 
assume t=0.4eV, U=7eV, and J=leV as typical values 
for transition metal oxidesjSi we obtain Ues/t=10 and 
J/t=2, where the system is expected to be near the phase 



boundary. Thus, the transition from gapful AFM to gap- 
less FM phases may be controlled by chemical pressure. 

In addition to the control of the Coulomb interactions 
or the bandwidth, it is possible that the magnetic field 
induces ferromagnetism in the Haldane systemi— The 
system is in the singlet ground state below a critical field 
due to the Haldane gap, indicating that magnetic mo- 
ment is suppressed. At the critical field, finite value of 
magnetic moment begins to grow. Note that in actual 
systems, spin-canted AFM phase occurs. Increasing fur- 
ther the magnetic field, the saturated FM state should 
occur, when we apply a high magnetic field corresponding 
to 4 Jafm"^ For instance, this value may be estimated as 
200 Teslas for LiVGe206, but smaller in other Haldane 
systems. The gapless excitation with respect to the or- 
bital is expected to be observed in a similar way that 
the orbital-wave excitation has been revealed by Raman 
scattering measurements in LaMnOa,^'' 

In summary, ground-state properties of the t2g-orbital 
degenerate Hubbard model have been investigated in or- 
der to clarify effects of orbital degeneracy in the spin 
S—1 Haldane system. In the AFM phase, the spin 5=1 
AFM Heisenberg model is reproduced, while in the FM 
phase, the system is described by the pseudo-spin T=l/2 
AFO Heisenberg model. Thus, we conclude that the low- 
energy physics is drastically changed due to the interplay 
of spin and orbital degrees of freedom. In particular, the 
gapless orbital excitation exists in the FM state. It may 
be interesting to observe such orbital excitation in the 
Haldane system under the high magnetic field. 
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